Critical fluctuations of the proton density in A+A collisions at 158^4 GeV 
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Studies of QCD suggest the existence of a critical point in the phase diagram of strongly interacting 
matter. Close to this point, according to recent theoretical investigations, the net-proton density 
carries the critical fluctuations of the chiral order parameter. Using intermittency analysis in the 
transverse momentum phase space of protons produced around midrapidity in the 12.5% most central 
C+C, Si+Si and Pb+Pb collisions at the maximum SPS energy of 158yl GeV we find evidence of 
power-law fluctuations for the Si+Si and Pb+Pb data. The fitted power-law exponent approaches 
the value expected for critical fluctuations. This suggests that the freeze-out states of these two 
systems are located in the phase diagram in the neighbourhood of the chiral critical point. 



PACS numbers: 25.75.-q 

The current understanding of the phase diagram of 
strongly interacting matter, as supported by several lat- 
tice QCD calculations [H, suggests the existence of a 
critical point (CP) at finite baryochemical potential and 
temperature. This CP is the endpoint of a line of first 
order transitions associated with the partial restoration 
of chiral symmetry when the temperature T, for given 
baryochemical potential \lb, increases beyond a critical 
value T c . For present experiments with heavy-ion col- 
lisions the search for the hypothetical "chiral" CP is of 



high priority. Many theoretical studies have been per- 
formed to find suitable observables 3 6]. To this end 
it is necessary to consider the fluctuations of the chiral 
condensate (qq) which is the order parameter of the chi- 
ral transition (g is the quark field). The quantum state 
which carries the quantum numbers as well as the criti- 
cal properties of the chiral condensate is the sigma-field 
a{x). In a heavy- ion collision the excitation of the vac- 
uum leading to the formation of the chiral condensate is 
likely to occur. However, the condensate will be unsta- 
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ble, decaying mainly into pions, at time scales character- 
istic of the strong interaction. The critical properties of 
the condensate are transferred to detectable tt + it~ pairs 
with invariant mass close to their production threshold 
5] . In a finite-density medium there is a mixing between 
the chiral condensate and the baryon density inducing 
critical fluctuations in the baryon density 0, H| as well. 
Furthermore, critical fluctuations of the chiral conden- 
sate are also directly transferred to the net-proton density 
through the coupling of the protons with the isospin zero 
CT-field [§j. Thus detecting the QCD CP through fluctu- 
ations of the net-proton density is a promising strategy. 

One of the main goal s of the NA49 experiment [l3| and 
its successor NA61 [TjJ is an extensive search for the CP 
by changing energy and size of the colliding nuclei in or- 
der to scan the phase diagram with a dense set of freeze- 
out states at different T and The observables of in- 
terest are either event-by-event fluctuations of integrated 
quantities like multiplicity fl2l ] and average transverse 
momentum [3] or local power-law fluctuations related to 
density-density correlations detectable through the mea- 
surement of factorial moments in momentum space 14 1 
within the framework of intermittency analysis 15[ . Our 
approach differs from the intermittency studies in the 
past [l6| in analysing quantities directly related to the 
order parameter of the phase transition (sigma conden- 
sate or net-baryon density). The advantage is that for 
a critical system the fluctuations of these quantities are 
expected to obey power laws with exponents determined 
by the corresponding critical exponents [4j, |6j . 

In the present work we calculate the second factorial 
moment of the proton density 
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as a function of the number M of subdivisions in each 
transverse momentum space direction (-1.5 GeV/c < 
Pt,x,Pt, v < 1.5 GeV/c), where rn is the number of pro- 
tons in the cell i. The analysed data were recorded by the 
NA49 experiment in A+A collisions at maximum CERN 
SPS energy of 158A GeV {^/sWn = 17.3 GeV). We can 
restrict the analysis to protons because the antiprotons 
are clearly far fewer for all considered systems. Theoret- 
ical work Q points out the emergence of intermittency - 
i.e. power-law behaviour of ^(M) as a function of M 2 - 
with expected exponent (intermittency index) </>2, C r = |j 
when the system freezes out exactly at the chiral CP. The 
critical intermittency index is determined from universal- 
ity class arguments. In the case of the two-dimensional 
system of interest (transverse momentum components) 
the above index can be revealed, if one studies protons in 
a small region (window) around midrapidity where their 
density reaches a plateau value. When moving the po- 
sition of the window towards beam rapidity, the critical 
behaviour is gradually masked [6|. Also, a weakening of 



the intermittent behaviour is expected as the distance of 
the freeze-out state from the CP increases, leading to a 
decrease of the value of^^and/or a modification of the 
power-law dependence [5J, [lj] . Our aim is to explore the 
systems C+C, Si+Si and Pb+Pb looking for power-law 
behaviour of the second factorial moment of the proton 
density in transverse momentum space as a signature for 
the approach to the CP. 

For the analysis we used data on the 12.5% most cen- 
tral collisions of "C", "Si" and Pb nuclei on C (2.4% in- 
teraction length), Si (4.4%) and Pb (1%) targets, respec- 
tively, collected by the NA49 experiment at ^Snn = 17.3 
GeV. The "C" beam as defined by the online trigger and 
offline selection was a mixture of ions with charge Z = 6 
and 7 (intensity ratio 69:31); the "Si" beam of ions with 
Z = 13,14 and 15 (intensity ratio 35:41:24) 0. The 
event statistics amounted to 210k events for "C" +C, 176k 
events for "Si"+Si and 1480k events for Pb+Pb. The 
standard event and track selection cuts of the NA49 ex- 
periment were applied as described in [l8j . Specifically, 
to avoid double counting split tracks we only accepted 
tracks for which the ratio of number of measured points 
to estimated maximum number of points in the TPCs 
exceeds 55%. Proton identification [l8| used the mea- 
surements of particle energy loss dE/dx in the gas of the 
time projection chambers. The inclusive dE/dx distribu- 
tion for positively charged particles in each reaction was 
fitted in 10 bands of momentum p to a sum of contribu- 
tions f a (dE/dx,p) from different particle species a with 
a = 7r, K, p, e. The probability P for a track with energy 
loss Xi and momentum p t of being a proton is then given 
by P = Pixi^/irixuPi) + f K (x uPi ) + fv(x hPi ) + 
f c (xi,pi)). The value of P for proton candidates had 
to exceed 80% for the "C"+C and "Si"+Si systems and 
90% for Pb+Pb. Moreover the candidates had to fulfill 
the criterion \dcm\ < Ubeam ~ 0.5 {ybeam ~ 2.9 in the 
cms) to avoid contamination by spectator protons [19|. 

In order to reveal a potentially underlying power-law 
behaviour of the factorial moments calculated from the 
data, we have to subtract the background due to uncor- 
rected protons (e.g. from hyperon decays) and/or the 
contamination by particles misidentified as protons. The 
background contribution is simulated by mixed events 
which, by construction, contain all the background effects 
but not the correlation signal. We then search for power- 
law behaviour of the corrected correlator Ai<2(M) = 



i<2 (M) — F2 m '(M) where the superscripts (d), (m) refer 



? {m) 



to data and mixed events respectively 

In Fig.[T]we show the second factorial moments F2(M) 
versus M 2 for the three analysed systems (red circles) . In 
the same figure we also plot the second factorial moments 
for the corresponding mixed events (black triangles). For 
the "Si" +Si and Pb+Pb systems we observe that for large 
M values the factorial moments of the data are clearly 
larger than those of the mixed events. The difference be- 
tween the two moments increases with increasing num- 
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ber of cells M 2 , a typical characteristic of intermittent 
behaviour. This observation is an indication for sizeable 
correlations among the produced protons. However, due 
to small statistics the errorbars of the F 2 values are quite 
large for "Si"+Si. In the "C"+C case the factorial mo- 
ment values of the data overlap with those of the mixed 
events suggesting the absence of an intermittency effect 
in this system. 
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FIG. 1: (Color online) Second factorial moments of the pro- 
ton density in transverse momentum space at midrapidity 
(—0.75 < ycM < 0.75) for the 12.5% most central collisions 
of (a) "C"+C, (b) "Si"+Si and (c) Pb+Pb at v /s7777 = 17.3 
GeV. The red circles (black triangles) represent the second 
factorial moment F^{M) of the data (mixed events) respec- 
tively. 

The background corrected correlator AF2 as a function 
of M 2 is shown in Fig. [5] The panels (a),(b),(c) illustrate, 
in log-linear scale, AF 2 (M 2 ) for "C"+C, "Si"+Si and 
Pb+Pb respectively. Only the region of large number of 
cells (M 2 > 2000) is displayed, where the intermittency 
effect, if present, is expected to show up [H, [2fjj . In 
order to make our analysis more quantitative we fitted 
the power- law function AF 2 (M) ~ (M 2 )* 2 shown by the 
lines in Fig. [2] A complication for the fitting procedure 
is the fact that the factorial moment values at consecu- 
tive M 2 are correlated. Several approaches like the corre- 
lated x 2 (Cx 2 ) [HI, the independent bin (IB) and the 
sparse binning (SB) [23| method were employed, which 
are designed to reduce the effect of these data point cor- 
relations. We found that for the Pb + Pb dataset the 
different methods give identical results fully compatible 
with the direct power-law fit (DPF). For "Si" + Si the 
results for (f> 2 have some spread (02, Cx 2 = 0.76 ± 0.28, 
fajB = 0.91±0.25 and 2 ,sb = 0.88±0.21), presumably 
due to small statistics. However, they are still compati- 
ble with each other within errors. The variation of these 
values is used to obtain a systematic error for "Si" + Si 
leading to: Sc/)^ = 0.12. Finally for the "C" + C sys- 



tem only the direct power-law fit is applicable. We also 
estimated the systematic error of 4>2 for the Pb + Pb sys- 
tem using different proton purity requirements (70%, 80 
% and 90%). Decreasing purity the power-law behaviour 
of the correlator is pushed to larger values of M 2 as ex- 
pected. Despite this, the associated power-law fits give 
(for the datascts with different purities) <j) 2 values which 
are compatible within errors confirming the efficiency in 
the subtraction of the background. Restricting the fit re- 
gion to the domain M 2 > 10000 we obtain (for the 90% 
purity case) <f) 2 = 0.36 ± 0.06 where we estimated the 
systematic error 54> 2 P ^y S = 0.06 as half the range of the 
three fit results. Due to low statistics it is not possible 
to increase purity for the "C" + C and "Si" + Si, but 
the combinatorial background is much smaller in these 
systems. 

For all performed fits the x 2 value per degree of free- 
dom was found in the range [0.4,0.9]. The quality with 
which the fit describes the variation of the correlator 
with the number of cells M 2 is measured in terms of 



the coefficient of determination R 2 24j . Using the DPF 
we obtained the intermittency index <f> 2 = —0.16 ± 0.44 
for the "C"+C system with R 2 « indicating that 
the "C" + C-correlator does not follow a power law in 
M 2 . The presence of many negative AF 2 values leads 
to a negative 4> 2 with a large error. For the "Si"+Si 
system the DPF gave cf> 2 = 0.88 ± 0.10 with R 2 = 0.60, 
well above the the minimum acceptable value of 0.5, 
supporting that the correlator in this case follows to a 
good approximation a power law. Finally for the Pb+Pb 
system the DPF gave cj> 2 = 0.37 ± 0.01 with R 2 = 0.96 
which is very close to the optimal value R 2 — 1 valid for 
an exact power-law dependence. 

Thus we observe significant power-law fluctuations of 
the proton density in transverse momentum space for the 
"Si"+Si and Pb+Pb systems at ^snn = 17.3 GeV with 
exponents approaching the QCD prediction (5/6). Both 
systems appear to freeze out near the CP in the phase 
diagram. It is illustrative to plot the fitted intermittency 
indices <f> 2 as a function of the size A of the nuclei (see 
Fig. [3]). To guide the eye we include in the plot also 
the QCD prediction for the CP 2 ,c<- = 5/6 (red solid 
line) [6j. Note that the 4> 2 value of the "C"+C system 
is added in the plot for completeness, although it is not 
associated with critical fluctuations since it is negative 
and the power-law behavior is questionable (R 2 « 0). 

An analogous intermittency effect was found recently 
[Til ] in central Si+Si collisions at ^snn = 17.3 GeV for 
7r + 7r~ pairs with invariant mass close to their production 
threshold (a- field fluctuations). The corresponding anal- 
ysis of Pb+Pb data was not possible due to limitations 
of the experimental resolution. 

Two remarks are in order before going to our conclu- 
sions. First, in Fig. U we show that the intermittency 
effect observed for the Pb+Pb system is confined, as the- 
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FIG. 2: (Color online) The correlator AF 2 (M) in log-linear 
scale, corresponding to the moments of Fig. 1, for the 12.5% 
most central collisions of (a) "C"+C, (b) "Si"+Si and (c) 
Pb+Pb at ^/sjviv = 17.3 GeV. The red solid lines are the 
results of DPF. 
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FIG. 3: (Color online) The dependence of the intermittency 
index <f>2 on the size A of the colliding nuclei. The red solid line 
indicates the QCD prediction at the CP: <j>2,cr = 5/6. Note 
that the plotted points were obtained by power-law fits which 
differ significantly in the value of the coefficient of determina- 
tion: R 2 C+C = 0.01, R 2 Si+S i = 0.60 and R 2 Pb+Pb = 0.96. 



oretically expected, to protons at midrapidity (the same 
is true also for "Si"+Si) where the density of the pro- 
duced fireball is approximately constant in rapidity @. 
We compared the correlator AF 2 for the Pb+Pb system 
calculated for protons with 1 < ycM < 1-4 (Fig. |4^) 
with that for -0.75 < y C M < 0.75 (Fig.gp). Clearly the 
effect gets suppressed for protons away from the midra- 
pidity region. We repeated the analysis for protons with 
-0.5 < y C M < 0.5 and found (j> 2 = 0.46 ± 0.04 with 
R 2 = 0.78. Thus within the errors there is no observ- 



able change of 4> 2 by further decreasing the size of the 
window around midrapidity. However, the R 2 value de- 
creases since the statistics becomes worse. 
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FIG. 4: (Color online) The correlator AF 2 for the 12.5% most 
central Pb+Pb collisions at ^/snn = 17.3 GeV using tracks 
with (a) 1 < ycM < 1.4 and (b) -0.75 < ycM < 0.75 plotted 
in a double logarithmic scale. 



Second, we calculated the correlator AF 2 for protons 
at midrapidity in the 12.5% most central Pb+Pb colli- 
sions at tJsnn — 8.8 GeV. In Fig. [5] we compare this 
result with that found at ^snn = 17.3 GeV. We clearly 
observe that the critical fluctuations disappear at lower 
energy. Thus, our analysis suggests that the critical bary- 
ochemical potential is closer to 240 MeV, the freeze-out 
baryochemical potential at ^/snn = 17.3 GeV, than to 
fi B ~ 380 MeV, the value at ^/sWn = 8.8 GeV [Hj]. 
A similar analysis cannot be performed for the "Si"+Si 
system at ^/snn — 8.8 GeV due to insufficient statistics. 

In summary, we performed a search for critical fluc- 
tuations in NA49 data using factorial moments of the 
proton density in transverse momentum space. We cal- 
culated the background subtracted second factorial mo- 
ment AF 2 in the 12.5% most central collisions of "C"+C, 
"Si"+Si and Pb+Pb at ^/sTHv = 17.3 GeV. We found 
a strong intermittency effect for the "Si"+Si and the 
Pb+Pb system at midrapidity. This effect is expressed 
through the power-law dependence: AF 2 (M) ~ (M 2 )^ 2 . 
The intermittency index <p 2 for the "Si"+Si system ap- 
proaches in size the QCD prediction for the CP. Our anal- 
ysis provides the first strong experimental indication for 
the existence of the QCD CP in the neighbourhood of the 
Si+Si and Pb+Pb freeze-out states in the phase diagram 
(T~ 150-160 MeV, fi B ^ 200-250 MeV). The absence 
of power-law correlations for Pb+Pb at ^/snn = 8.8 GeV 
(fiB ~ 380 MeV) excludes the location of the CP close 
to the value of fis predicted by the lattice calculations 
of Ref.[2j. The exact location of the CP may possibly 
be further constrained by measurements of proton den- 
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FIG. 5: (Color online) Correlator AF 2 for the 12.5% most 
central Pb+Pb collisions at ^snn = 17.3 GeV (red circles) 
and ^/sjvjv = 8.8 GeV (blue triangles) in the midrapidity 
region (-0.75 < ycM < 0.75). 

sity fluctuations in the experiment NA61 studying A+A 
collisions with small and intermediate size nuclei. Signif- 
icant new results on fluctuations may also be expected 
from the RHIC Beam Energy Scan program [26j . 
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